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1.0 Summary 

An ultrasonic scan procedure using the pulse-echo contact 
configuration was employed to obtain maps of pore fraction 
variations in sintered silicon nitride samples in terms of ultrasonic 
material properties. Ultrasonic velocity, attenuation coefficient, 
and reflection coefficient images were obtained simultaneously 
over a broad band of frequencies (e.g., 30 to 1 10 MHz) by using 
spectroscopic analysis. Liquid and membrane (dry) coupling 
techniques and longitudinal and shear-wave energies were used. 
The major results include the following: Ultrasonic velocity 
(longitudinal and shear wave) images revealed and correlated 
with the extent of average through-thickness pore fraction varia- 
tions in the silicon nitride disks. Attenuation coefficient images 
revealed pore fraction nonuniformity due to the scattering that 
occurred at boundaries between regions of high and low pore 
fraction. Velocity and attenuation coefficient images were each 
nearly identical for machined and polished disks, making the 
method readily applicable to machined materials. Velocity images 
were similar for wet and membrane coupling. Maps of apparent 
Poisson’s ratio constructed from longitudinal and shear-wave 
velocities quantified Poisson’s ratio variations across a silicon 
nitride disk. Thermal wave images of a disk indicated transient 
thermal behavior variations that correlated with observed varia- 
tions in pore fraction and velocity and attenuation coefficients. 

2.0 Introduction 

Traditional nondestructive evaluation (NDE) imaging meth- 
ods applicable to the characterization of high-performance 
materials, such as vehicular engine ceramics and composites, 
include x-radiography, ultrasonic imaging, and thermal inspec- 
tion (ref. 1 ). In their standard configurations, these methods 
measure an amplitude of transmission of the respective energy 
used. These methods provide an efficient, quasi-quantitative 
measure of material homogeneity but may lack the precision 
necessary for microstructure evaluation for high-performance 
materials (refs. 2 and 3). The use of materials in demanding 
high-performance applications will require detailed quantitative 
knowledge of microstructural and compositional variability for 
defining acceptable levels of variability and for rejecting those 
materials and processes that yield sample-to-sample and within- 
sample variations likely to result in unacceptable property (e.g., 
strength, thermal conductivity, or oxidation resistance) varia- 


tions. Variability must be precisely characterized either directly 
in terms of property measurement or indirectly through micro- 
structural characterization where microstructure-property rela- 
tions have been previously established. In support of the latter, 
quantitative NDE methods for the determination of the global 
microstructural state are currently being developed and refined 
to aid in advanced materials development efforts (refs. 4 to 7). 

Repeated, uniformly spaced ultrasonic contact measurements 
have been successful in previous studies for quantifying and 
mapping inhomogeneity in various ceramics (SiC, AI 2 O 3 , 
YBa 2 Cu 307 ) and metals in terms of ultrasonic material proper- 
ties (refs. 6 and 8 to 12). Such ultrasonic material properties 
include velocity, reflection coefficient, and attenuation coeffi- 
cient. In most cases, velocity mapping has been the preferred 
ultrasonic property map due to its repeatable nature, relatively 
low sensitivity to material surface condition and edge effects 
(ref. 13), and relative ease of quantitatively interpreting results in 
terms of microstructural features (ref. 14). Spatial microstructural 
variations previously characterized in terms of ultrasonic velocity 
include variations in pore fraction, granular texture, and residual 
stress (refs. 9, 10, and 12). Reflection coefficient and attenuation 
coefficient maps derived from contact scanning have been used 
to a lesser extent for characterizing spatial variations in polycrys- 
talline materials (ref. 6 ) because of an expected sensitivity to 
surface condition (ref. 15) and edge effects (ref. 13) and more 
difficulty in quantitatively interpreting results in terms of 
microstructrual features (ref. 6 ). 

In this study, the ultrasonic contact pulse-echo method was 
used in an automated scanning configuration in an attempt to 
quantify and map nonuniformity in polished and machined 
sintered silicon nitride (S^N^ disks containing distinct regions of 
high pore fraction formed during sintering. Here, the configura- 
tion was significantly (hardware and software) modified from 
that of previous approaches (refs. 9, 10, and 16). Modification 
allowed the examination of any two-dimensional, complex- 
shaped (e.g., circular, tensile shaped) sample by using small- 
diameter probes. The modification also allowed the use of dry as 
well as wet coupling and thus shear as well as longitudinal wave 
energies. Reflection coefficient, attenuation coefficient, and 
velocity maps of the S^N 4 disks were formed over a wide 
frequency range (e.g., 30 to 110 MHz) by using spectral analysis 
after acquisition of ultrasonic waveform data. These three types 
of ultrasonic property maps were examined to determine which 
ultrasonic property map is most effective in characterizing 
nonununiformity and the nature of the correlation, if any, between 
different ultrasonic property maps. Spectral analysis was employed 
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to determine if inhomogeneity was a function of frequency and to 
characterize elastic (dispersion) behavior for Additionally, 

maps of Poisson’s ratio for a disk were constructed from longitu- 
dinal and shear-wave velocity images. Thermal wave imaging 
and Vickers microhardness measurements were used to deter- 
mine if physical behavior variations correlated with pore fraction 
variations that had been characterized metallographically and in 
ultrasonic images. A positive correlation would indicate the 
usefulness of the ultrasonic measurement as an aid for predicting 
physical behavior. 

The authors wish to thank the Metallographic Support Staff for 
their help in sample preparation and analysis, James R. Bodis for 
his help in obtaining ultrasonic C-scan and microfocus radio- 
graphy images, and Maurice Bales of Bales Scientific, Inc., and 
Elliot Cramer of NASA Langley Research Center for their help 
in obtaining thermal wave images. 


3.0 Background 

3.1 Importance of S 13 N 4 

The past 20 years have seen significant interest in ceramic 
materials for high-temperature structural applications, such as 
vehicular engines (refs. 17 and 18). Si 3 N 4 is a candidate material 
for such applications because of its low density, high melting 
point, high strength, and resistance to wear, thermal shock, and 
creep (refs. 17 and 18). Si 3 N 4 has application as a monolithic 
material, a toughened monolithic material, and the matrix mate- 
rial in ceramic matrix composites (refs. 17 to 19). The Si 3 N 4 used 
in this study can also be considered a model material for demon- 
strating the capabilities of the modified ultrasonic contact scan 
technique. 



(a) Diagram of buffer rod-couplant-sample pulse-echo contact 
configuration. 



(b) Resulting waveforms for pulse-echo contact technique. 


3.2 Relationships Between Ultrasonic Properties and 
Experimental Variables 

The significant material-behavior-related variables in this study 
were ultrasonic frequency and pore fraction. A detailed discus- 
sion follows concerning expected relationships between these 
variables and the ultrasonic reflection coefficient, attenuation 
coefficient, and velocity for Si 3 N 4 when the contact pulse-echo 
experimental technique (refs. 20 and 2 1 ) is used. In this technique, 
an ultrasonic pulse generated by a piezoelectric crystal is propa- 
gated into a buffer rod-couplant-sample configuration (figs. 1(a) 
and (b)). In this investigation, the buffer rod, the couplant, and the 
sample are fused silica (SiC> 2 ), diffusion pump fluid, and Si 3 N 4 , 
respectively. The ultrasound reflects off the front and back 
surfaces of the sample and travels back to the crystal. Figures in 
the immediate upcoming sections show the expected frequency 
and pore fraction dependence of velocity, reflection coefficient, 
and attenuation coefficient for the following conditions: contact 
pulse-echo configuration, frequency range of 30 to 110 MHz, 
pore fraction range of 0 to 0.1, poreradius of 1 fim , and ultrasonic 
velocity in fully dense Si 3 N 4 of 1.12 cm/jLisec. 



(c) Schematic (top view) of ultrasonic contact scan procedure 
showing examples of successive transducer positions along 
X- and /-dimensions of sample. 

Figure 1 . — Ultrasonic measurement method (FS = front-surface 
reflection; B 1 = first back-surface reflection; B 2 = second 
back-surface reflection). 
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Impinging 

ultrasonic 

Schematic of volume of pulse Thickness 



(d) Schematic (three-dimensional view) showing volume of 
sample ultrasonically interrogated at each scan point for 
velocity and attenuation coefficient measurements. 

Figure 1 . — Concluded. 

Velocity. — Brittle materials, such as monolithic ceramics, 
generally exhibit linear elastic behavior at room temperature 
(ref. 22), and therefore dispersion (velocity varying with fre- 
quency, ref. 23) and frequency-dependent attenuation are not 
expected to be present to the extent that they would be for 
viscoelastic media, such as polymeric materials, where energy 
absorption is the dominant attenuation mechanism (refs. 24 and 
25). However, the presence of defects in ceramics, such as pores, 
grain boundaries, and dislocations, can lead to frequency- 
dependent effects because ultrasonic energy is scattered from 
such defects. For the case of multiple scattering where the 
scatterers are cavities, a low cavity concentration exists, and 
Rayleigh scattering predominates, Sayers has shown the follow- 
ing relationship for velocity (refs. 26 and 29): 

V= V o {l + 0.5P[A + B(*a) 2 ]} (1) 


frequency by using equation (2). The velocity Vsi 3 N 4 is expected 
to linearly decrease with increasing pore fraction (from 
1.12 cm/fuscc at full density to about 1.06 cm/|isec for 0.1 pore 
fraction) for all frequencies and to show little dependence on 
frequency over the specified frequency range. Figure 2(b) shows 



(a) Theoretical velocity versus pore fraction for various 
frequencies (eq. (2)). 



Pore fraction, P 

(b) Measured velocity versus pore fraction where ^SioNy* " 
-0.01 4 x Percent porosity + 1.12 was previously 3 4 
observed empirically (0.98 correlation coefficient; ref. 14). 


where A and B are constants that depend on the ratio of longitu- 
dinal to shear-wave velocity, ka-lnaflV, P is pore fraction, V Q 
is the velocity in a fully dense material containing no pores, V is 
the velocity in a material containing pores, a is pore radius, and / 
is frequency. For the Si 3 N 4 in this study, VJVs ~ 1 .8; thus^4 = -1.1 
andZ? = — 1.8 according to Sayers’ analysis (ref. 26). Equation (1) 
thus becomes 


v Si,N, =U2 cm //i sec 


1 + 0.5P 


— 1.1 — 1.8 


2 nfa 


1.12 cm / fi sec 


( 2 ) 


Figures 2(a) and (c) show the expected theoretical dependence 
of the ultrasonic velocity in Si 3 N 4 (Vsi 3 n 4 ) on pore fraction and 


1.12 

1.10 

1.08 
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(c) Theoretical velocity versus frequency for various pore 
fractions (eq. (2)). 


Figure 2. — Expected dependence of ultrasonic velocity in Si 3 N 4 
on pore fraction and frequency. 
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velocity versus pore fraction for where the relationship 

Vsi 3 N 4 = ~ 0.014 x Percent porosity + 1.1 2 was previously ob- 
served empirically (0.98 correlation coefficient) in reference 14. 

Reflection coefficient. — The theoretical expression for the 
buffer rod-couplant-sample (BCS) pressure reflection coeffi- 
cient is (ref. 28) 


( 7 ^ 

2 3 J 

COS 

~2nf{8X) 

v 2 

+ j 

V 

k Z 3 Z 2 j 

sin 

2nf(8X)~ 

V 2 

f 7 ^ 

1+4. 

Z 3 ) 
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'2 nf(SX) 
y 2 

+ j 

% + V 

< Z 3 Z 2 ) 

sin 

'2nf{SX) 

. V 2 . 


(3) 

where / is the frequency of the ultrasonic transducer; V 2 is the 
velocity of sound in the couplant layer; Z\, Z 2 , and Z 3 are the 
acoustic impedance of the buffer rod, couplant, and sample, 
respectively; and fiYis the couplant thickness. The magnitude of 
this complex expression is (ref. 15) 


l*l= 



~2jtf(8X) 

. v 2 . 

<n| m 

N 1 N 
+ 

7 ^ 
L \ 

Z 2 y 
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. 2 
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'2 */(»)" 
. V 2 . 


( 7 V 

1 H COS 2 

A z 3 ) 

. V 2 . 

Ub 

7 ^ 

+ 4. 

z 2 J 

2 

. 2 

sin 

’2 nf(SK) 

. V 2 . 



(4) 

The acoustic impedance of the sample Z 3 is given by (ref. 13) 

Z 3 =p 3 V 3 (5) 

where P 3 and V 3 are the sample density and velocity, respectively. 
As the number of scatterers near the material surface increases, it 
is expected that P 3 and V 3 will be altered, changing Z 3 and thus 
altering R. The manner in which R is expected to change with 
material change depends on the compositional nature of the 
scatterer (pore, second-phase inclusion, etc.). For example, if the 
scatterers are pores, Z 3 is expected to decrease with increasing 
volume fraction of pores as P 3 and V 3 are effectively decreased 
according to (ref. 14) 

2 3 =P0 >3 (1 -*)%(!-*) (6) 

where P 03 , Vo, 3 , and P are the density of a sample with no pores, 
the velocity in a sample with no pores, and the pore fraction in 
the sample, respectively. Substituting the right side of equa- 
tion ( 6 ) into equation (4) for Z 3 and substituting the values 
SX = 1 pm, Z S jo 2 =l-31 g/cm-sec, Z coup | ant = 0.15 g/cm-sec, 


Po,Si,N 4 =3.3 g/cm 3 , and V 0 > si 3 N 4 = 1.12 g/cm 3 in equation (4) 
allows the magnitude of the reflection coefficient /?si 3 N 4 for the 
buffer (SiC> 2 ) — couplant (Dow Coming 704) -sample (S 13 N 4 ) 
configuration used in this study to be expressed in terms of 
frequency and pore fraction. Figures 3(a) and (b) show the 
expected dependence of the theoretical reflection coefficient 
/?Si^N 4 on P ore fraction and frequency for the BCS configuration 
of this study (couplant thickness, 1 pm). The reflection coeffi- 
cient /?si 3 N 4 is expected to linearly decrease with increasing pore 
fraction for all frequencies and more steeply at lower frequen- 
cies and to increase with increasing frequency for all pore 
fractions. 

Changing the pore fraction may also affect the average surface 
roughness, which may then alter the average couplant thickness 
bX and change ^si 3 N 4 - Figure 3(c) shows the frequency depen- 
dence of /?si 3 N 4 f° r several different couplant thicknesses for 
frequencies from 30 to 110 MHz. For a given frequency, /?$i 3 N 4 
increases with increasing couplant thickness. For a given couplant 
thickness, /?si 3 N 4 increases with increasing frequency except for 
larger couplant thicknesses such as 5 pm, where /?si 3 N 4 initially 
increases slightly with increasing frequency and then begins to 
decrease slightly with increasing frequency at about 75 MHz. 

Attenuation coefficient. — A theory has been developed for 
attenuation due to grain scattering from equiaxed grains in a 
single-phase, fully dense polycrystalline solid. This theory pre- 
dicts that attenuation is a function of average grain size D , sonic 
wavelength A, and frequency / and that this functionality depends 
on the ratio of sonic wavelength to grain size. Table 1 reviews the 
frequency dependence of the attenuation coefficient and the 
scattering mechanisms for polycrystalline materials (ref. 13). 
Most materials have a range of grain sizes and contain additional 
scatterers besides grain boundaries, such as pores, second-phase 
inclusions, and dislocations, that affect the attenuation coefficient. 

The two-phase scattering theory developed by Truell for the 
case where scattering from individual scatterers is independent 
relates the frequency and the size and number of spherical 
scatterers to the attenuation coefficient a according to (refs. 27 
and 29) 

cc = ^n oY (7) 

where n 0 is the number of scatterers per unit volume and y is the 
scattering cross section, which is a function of frequency and 
scatterer radius (refs. 27 and 29). The normalized scattering cross 
section yv is defined as (refs. 27 and 29) 

Yn= JL J ( 8 > 

na 

where a here is the scatterer radius. Reference 27 used the two- 
phase scattering theory to show that the attenuation coefficient 
due to scattering from cavities and inclusions increases with 
increasing frequency. 
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Theoretical reflection coefficient, R 



(a) Reflection coefficient versus pore fraction for various 
frequencies. (Couplant thickness, 1 pm.) 



(b) Reflection coefficient versus frequency for various 
pore fractions. (Couplant thickness, 1 pm.) 



(c) Reflection coefficient versus frequency for various 
couplant thicknesses (eq. (4)). 


Figure 3. — Expected dependence of theoretical BCS reflection 
coefficient on pore fraction, frequency, and couplant thickness 
for buffer of Si0 2 , couplant of Dow Coming 704 oil, and Si 3 N 4 
sample. 


TABLE 1.- FREQUENCY DEPENDENCE OF 
ATTENUATION COEFFICIENT AND 
SCATTERING MECHANISMS FOR 
POLYCRYSTALLINE MATERIALS 
[Experimental constants C, are dependent on grain, 
geometric and elastic anisotropy, density, and 
longitudinal and shear-wave velocity factors.] 


Wavelength/ 
grain size 
relation 

Attenuation 

mechanism 

Attenuation 

coefficient 

Independent 
A »D 
k~D 
X«D 

True absorption 
Rayleigh scattering 
Phase scattering 
Diffusion scattering 

O-a-CJ 
a R = C R DY 
a p = CpDf 1 
a d -C d !D 


Reference 30 has calculated relative changes in the normalized 
scattering cross section due to scattering from pores in Si 3 N 4 as 
a function of wavenumber times scatterer radius ka, as shown in 
figure 4(a). Since ka = (2/ra/)/(Vsi 3 N 4 )> where Vsln 4 is the 
velocity of ultrasound in relative change in the normalized 
scattering cross section is a function of frequency for a given pore 
radius and Vsj 3 N 4 - Thus, it is possible to estimate changes in the 
Si 3 N 4 attenuation coefficient as a function of pore concentration 
(number of pores per cubic centimeter) and changes in frequency 
by using the following methodology and figure 4(a): 

(1) Assume a constant Vsi 3 N 4 of 1.08 cm/jisec and an average 
pore radius a of 1 fim (or 0.0001 cm). Then, the x axis in 
figure 4(a) becomes (5.8X 10 -4 sec)/, where/is in megahertz. 

(2) Select a frequency range over which it is desired to deter- 
mine attenuation coefficient changes (e.g., 30 to 110 MHz). 

(3) From the curve shown in figure 4(a), relative changes in 
normalized scattering cross section can be determined for changes 
(increases) in frequency, and a relationship can be derived between 
A yn and A/. A relationship was derived by using nonlinear 
regression analysis with a power function model 

A y N = ^A/^j + Z^, forcing the £>2 exponent to be 0.5 (the 

curve shape in fig. 4(a) approximates y^=(ka)°‘ 5 ). For an / of 
40 to 110 MHz (corresponding to a A / of 10 to 80 MHz from a 
reference frequency of 30 MHz), 

A y Ns ^ t =0.3396(A/ 05 )- 0.6352 (9) 

where A y^ is in nepers ( 1 neper = 8.686 dB) and A/is in megahertz. 
This relationship, shown in figure 4(b), has a coefficient 
of determination of 0.9996 and a final sum of squares residuals 
of 0.0013. 

From equation (7), for constant n () , changes in attenuation 
coefficient as a function of changes in scattering cross section can 
be expressed as 

= 0 °) 


5 



(a) Relative changes in normalized scattering cross section 
due to scattering from pores in Si 3 N 4 as a function of 
wavenumber times scatterer radius ka. From reference 30. 



40 50 60 70 80 90 100 110 

Actual frequency, f, MHz 

(b) Expected changes in normalized scattering cross section as 
function of change in frequency from 30 MHz (corresponding to 
frequency of 40 to 1 1 0 MHz for pore radius of 1 pm and velocity 
in Si 3 N 4 ( V Sj N ) of 1 .08 cm/psec (derived from fig. 4(a))). 


Rearranging equation (8) and substituting for yin equation (10) 
gives 


Aa = 


10 8 cm 2 
2 


nn Q Ay N 


( 11 ) 


Substituting the right side of equation (9) into equation (11) gives 


A «Si 3 N 4 = 


^10 8 cm 2 


;r« o [(0.3396)A/°- 5 -0.6352] (12) 


where Aasi 3 N 4 * s in nepers per centimeter and A/is in megahertz. 
This relationship ignores the scattering conditions and frequency 
relationships shown in table I and just takes into account the shape 
of the curve in figure 4(b). From the ratio of the sonic wavelength 
(-200 fj m) to the average scatterer (grain and pore) size (-1 /jm) 
observed for the Si 3 N 4 used in this study, one would predict 
Rayleigh scattering ( a proportional to f 4 ). Figure 4(c) shows 
the expected dependence of changes in a$i 3 N 4 (A«si 3 N 4 ) on 
positive changes in frequency A/ from 30 MHz, assuming that 
pore scattering is the sole cause of attenuation for several different 
values of pore concentration (number of pores per cubic centime- 
ter) n Q . The «si 3 N 4 is expected to increase with increasing fre- 
quency for any pore concentration, but the change only becomes 
significant for n 0 > 10 8 pores/cm 3 . 

For nonindependent scattering (i.e., multiple scattering condi- 
tions where the scatterers become sufficiently dense and closely 
spaced and where the wavelength becomes comparable to the 
scatterer size), equation (7) cannot be used to predict the attenu- 
ation coefficient (ref. 29). 



1 1 1 1 I I 1 1 

40 50 60 70 80 90 100 110 

Actual frequency, f, MHz 

(c) Expected dependence of changes in a Si m ( Aa Si 3 N 4 ) 
due to scattering from pores on positive changes in 
frequency from 30 MHz for several pore concentrations 
(eq. (12)). 

Figure 4.— Expected changes in normalized scattering cross 
section and attenuation coefficient for Si 3 N 4 . 


3.3 Ultrasonic Measurements 

Velocity. — Velocity can be determined in several ways by using 
the pulse-echo configuration. In this study, cross-correlation and 
phase methods were used by applying frequency-domain signal 
processing techniques. In the cross-correlation method, the time 
delay ^between the back-surface-reflected pulses# \(t) and #2(0 
(figs. 1(a) and (b)) was obtained by cross-correlating the pulses 
with digital signal processing (ref. 21), and the velocity V was 
obtained from r^and the sample thickness X according to (ref. 21) 

V = 2— (13) 

^ O 

Velocity determined from cross-correlation is essentially a 
pulse velocity because the entire wave train (containing a broad 
band of frequencies) is considered in the calculation. Cross- 
correlation is a preferred method for determining velocity, since 
it produces accurate velocities even with noisy signals (ref. 21). 
If pulse smearing occurs owing to dispersion- and/or frequency- 
dependent attenuation, phase methods (ref. 3 1 ) where the velocity 
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is determined for a particular frequency may yield more accurate 
results. The frequency-dependent phase velocity is given by 
(ref. 13) 

= , A0=0,-0 2 (14) 

A 6 

where 


0j(/)=tan 1 
and 


0 2 (/)=tan 1 


where /is frequency and B\(f) and Z?2(/) are the Fourier transfor- 
mations of the time-domain pulses B\(t) and Z?2(0> respectively. 

Attenuation coefficient. — The frequency-dependent attenua- 
tion coefficient for the pulse-echo configuration is given by 
(refs. 6 and 13) 


Re[B 2 (/)]j 


(16) 



a(f) = — In 
v ' 2X 


h(/)IMI 


(17) 


where | R(f) | is the frequency-dependent reflection coefficient 



(18) 


/is frequency, I FS](f) | and \FS2(f)\ are the Fourier magnitude 
spectra of the time-domain pulses FS\(t) and FS2(t) without and 
with the sample present on the buffer rod, respectively, and 
I B\(f) | and 1 82 (f) I are the Fourier magnitude spectra of the 
time-domain pulses B\(t) and #2(0> respectively. 

The errors associated with velocity, reflection coefficient, and 
attenuation coefficient measurements obtained by using the con- 
tact pulse-echo method are discussed in appendix A. 


4.0 Experiment 

4.1 Si 3 N 4 Samples 

In this investigation, seven disk-shaped samples (labeled 

SN 1 to SN7) approximately 3.5 cm in diameter and 0.3 cm thick 
were manufactured. Details of the fabrication process are given in 
appendix B. The Si^ISL* disks were 3.24 ± 0.03 g/cm 3 (> 98 per- 
cent theoretical density) as determined from immersion density 


measurements. From quantitative image analysis at 400 X and 
800 X, pore dimensions were found to range from <1 /im (~90 per- 
cent of the pores) to as large as 10 fim with an average of 
approximately 1 |tm. Grain dimensions for representative plasma- 
etched (ref. 32) sections of samples ranged from approximately 
0.5 to 10 //m, with an average grain dimension of approximately 
ljUm from scanning electron microscopy (SEM) analysis at 
2000 X to 10 000X. 

At the locations where disk separators (setters) had touched the 
disks during sintering, visibly lighter colored regions (relative to 
the surrounding regions on the disks) were noted. In some cases, 
three such regions were apparent on each side of the disk (fig. 5). 
For one of the disks, no such regions were visible. It was deter- 
mined that these lighter areas contained a significantly higher 
volume fraction of pores than did the surrounding regions (fig. 6). 
Also, the central sections of the disks were generally more porous 
than the edge regions. Figures 7(a) and (b) show the percent of 
porosity on the surface and the average percent of porosity for 
lines through a section thickness, respectively, from quantitative 
image analysis at 800 X of disk SN2 for various locations. 
Figures 8(a) and (b) show representative profiles of percent of 
porosity versus depth through the thickness for the lighter colored 
and darker colored surrounding regions, respectively. 

Energy-dispersive and wavelength-dispersive x-ray scans 
showed no chemical (silicon, nitrogen, oxygen, yttrium, or boron) 
differences between the high-pore-fraction region and the sur- 
rounding regions. SEM analysis revealed no obvious grain size 
distribution differences between the high-pore fraction region 
and the surrounding regions. For one sample area examined in 
disk SN2, an yttria-rich glass phase was identified at grain 



Figure 5. — Disk SN6 showing three visibly lighter colored regions 
(relative to surrounding regions). The other side of this disk 
showed the same pattern but offset 60°. 
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(a) Pore fraction on disk surface. 



(b) Surrounding darker colored, less porous region. 
Figure 6.— Photomicrographs of disk SN2. 
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(b) Average through-thickness pore fraction at selected 
positions of disk SN2. 


Figure 7. — Pore fraction profiles from quantitative image 
analysis at 800x of disk SN2 for various locations. 


boundaries that was approximtely 30 percent less prevalent in the 
high-pore-fraction region (-0.09 volume fraction glass-rich phase) 
than in the surrounding region (-0.12 volume fraction glass-rich 
phase). For the same area, an approximately 5 X difference in pore 
fraction was noted between the high-pore-fraction region 
(-0.05 pore fraction) and the surrounding region (-0.01 pore 
fraction). In some cases, a glass-like material was observed at the 
bottom of and lining the pores. This material was not identified by 
using x-ray diffraction and energy -dispersive spectroscopy. A 
trace amount of an iron silicide phase was also identified at the 
grain boundaries. Figure 9 shows typical surface profiles for 
machined and polished surfaces of the Si 3 N 4 disks. Machined 
disk surfaces were characterized by less than 1-jUm peak-to-peak 
roughness. No differences in surface condition were observed 
between polished regions at a disk center and at a high-pore- 
fraction area. 


42 Ultrasonic Contact Scanning 

Ultrasonic measurements using the pulse-echo contact tech- 
nique were obtained over an ordered array of points across the 
entire surface of each S^N 4 sample disk by means of an auto- 
mated scanning technique (fig. 1). A single transducer in com- 
bination with spectral analysis allowed complete material 
characterization over a wide band of frequencies. This method 
results in considerable time savings over scanning using discrete 
frequency techniques. The scan technique was hardware and 
software modified from that of previous approaches (refs. 9, 10, 
and 16) to allow the scanning of entire samples without probe 
damage due to movement over sample edges (ref. 16), the use of 
transducer probes having buffer rod faces as small as 0.25 in. 
(0.6 cm) in diameter, and the use of membrane (dry) coupled 
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(b) Darker colored, less porous surrounding region. 


Figure 8. — Pore fraction profiles from quantitative image analy- 
sis at 800x of pore fraction versus depth through thickness for 
disk SN2. 


(allowing shear-wave scans) as well as liquid (wet) coupled 
measurements. 

The sample disks were mounted in a Lucite holder in prepara- 
tion for scanning. Each scan consisted of a 47- by 41 -grid of 
measurements with each measurement separated from the next by 
1 mm. Scans were run twice in most cases to verify reproducibil- 
ity. Scanning was accomplished through the use of computer- 
controlled x-, y -, and z-micropositioner tables that were used 

(1) to position the sample and holder under the transducer probe, 

(2) to raise the sample and holder so that the probe was in contact 
with the sample (the contact force on the probe for each measure- 
ment in a scan was 1 2 ± 0. 1 lb for wet-coupled scans and 20 ± 0. 1 lb 
for dry -coupled scans; the couplant layer thus should be approxi- 
mately uniform for all measurements in the scan), (3) to bring the 
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(c) Polished, at lighter colored, more porous region. 
Figure 9. — Typical surface profiles for surfaces of Si 3 N 4 disks. 


sample away from the probe after the measurement was made, 
and (4) to reposition the sample in the x- or y-direction (1-mm 
increment) for the next measurement. Echo waveforms from the 
sample front and back surfaces were digitized (512-point hori- 
zontal and 5 1 2-point vertical resolution), averaged over 64 wave- 
forms, and stored. Each waveform was digitized at the voltage 
gain setting that allowed use of the full vertical resolution capa- 
bility of the digitizer. 

Digitized waveforms were Fourier transformed and subse- 
quently processed to yield ultrasonic properties. Beam-spreading 
diffraction corrections were made to the amplitude spectra of each 
echo according to the method of Papadakis (ref. 33). (The effect 
of diffraction on the phase of the ultrasonic waveforms is small, 
and no corrections to phase were made (ref. 24). Ultrasonic 
images were constructed over a wide band of frequencies (e.g., 30 
to 110 MHz, at 10-MHz intervals) from reflection coefficent, 
attenuation coefficient, and phase velocity values obtained at 
each scan point. One cross-correlation velocity image for each 
scan was also constructed. The velocity and attenuation coeffi- 
cient images can be thought of as two-dimensional projections 
representing averaged microstructural information for the vol- 
ume of sample scanned (fig. 1(d)). A computer program was 
written to allow continuous image formation at discrete frequen- 
cies and subsequent storage on an optical disk. Acontinuous scale 
consisting of 256 shades of gray (or color) and linear interpolation 
between points allowed the display of spatial variations in these 
ultrasonic properties. A data clipping or filtering scheme was 
performed to remove faulty property values resulting from incor- 
rect digitization and edge effects. Failure to apply the filtering 
resulted in misleading images and poor contrast between non- 
uniformities and surrounding regions in the disk. 

Transducers used were broadband 10-MHz longitudinal wave, 
5-MHz shear wave, 50-MHz longitudinal wave, and 125-MHz 
longitudinal wave with the megahertz value denoting the nominal 
center frequency. The 10-MHz longitudinal and 5-MHz shear- 
wave transducers consisted of a critically damped lead metaniobate 


9 


crystal with an active diameter of 0.6 cm bonded to a polystyrene 
buffer rod 0.6 cm in diameter and 0.8 cm in length. A thin 
(0.04 cm) pliable polymer material attached to the buffer rod 
allowed ultrasound to be transmitted into the sample without the 
need for liquid coupling (Ultran, Boalsburg, PA). The 50- and 
125-MHz longitudinal wave transducers consisted of a lead 
metaniobate crystal with an active diameter of 0.3 cm bonded to 
a silica glass buffer rod 0.6 cm in diameter and 0.8 cm in length 
(Panametrics, Waltham, MA). Liquid couplant (diffusion pump 
fluid, 704 type, Dow Coming, Midland, MI) was necessary 
between the transducer probe and the sample to allow ultrasonic 
transmission. Contamination from the couplant was assumed to 
be minimal because of the high average density (>98 percent) of the 
samples. Pulser receivers of 0.01 to 10 MHz and 1- to 150-MHz 
bandwidth were used with the 5- and 10-MHz and 50- and 125- 
MHz transducers, respectively. Signal-to-noise ratios were high 
(-250) for the wet-coupled configuration and moderate to low 
(-10 to 50) for the dry -coupled configuration. 

43 Further Characterization 

Ultrasonic immersion C-scans. — Immersion (water) C-scans 
were obtained by using 10- and 35-MHz focused transducers to 
qualitatively examine the uniformity of the Si 3 N 4 disk samples 
and to compare them with contact scan images of the samples. 
Data were acquired every 100 jum in both the jc- andy-directions. 
Transmission amplitude images were obtained by gating front 
and back surface reflections, and transmitted pulses were obtained 
by using pulse-echo and through-transmission techniques, 
respectively. Scans were run with the polished faces facing the 
sending transducer. The images show the amplitude of ultrasonic 
transmission in terms of gray scale with lighter gray scale corre- 
sponding to higher ultrasonic transmission. Images were contrast 
expanded for optimum viewing of nonuniformity. 

X-ray radiography- M icrofocus rad iography was used to qual i- 
tatively examine the uniformity of the disk samples and to 
compare them with contact scan images of the samples. Radio- 
graphs were made through the thickness of the disks with conven- 
tional contact methods under conditions optimized to yield high 
accuracy and contrast. The source-to-disk distance was 
approximately 18 in., and the exposure conditions were 35 kV, 
10 mA, and 4 min. The images show x-ray attenuation in terms of 
gray scale, with lighter gray scale corresponding to higher 
attenuation. 

Thermal wave imaging- Thermal wave imaging was used to 
determine if thermal behavior variations correlated with pore 
fraction variations that had been characterized metallographi- 
cally and in ultrasonic images of disk SN6. The polished surface 
of disk SN6 was coated with a high-emissivity material to allow 
maximum absorption of thermal radiation into the sample. Four 
flashlamps were located at 45° angles to, and 18 in. away from, 
the polished surface. The lamps were used to externally heat SN6. 
Each flashlamp’s power was 1.5 kJ. A high-resolution infrared 
camera was located in front of the disk on the same side as the flash 
lamps, and data were obtained at a rate of 20 frames per second. 


The images show thermal energy density (temperature) in terms 
of gray scale, with lighter gray scale corresponding to higher 
thermal energy density (temperature). 

Microhardness measurements .-Vickers microhardness meas- 
urements were used to determine if hardness variations correlated 
with port fraction variations that had been characterized metallo- 
graphically and in ultrasonic imagesofdisk SN3. Microhardness 
measurements were performed at the regions of high pore fraction 
and at surrounding regions of average pore fraction in a polished 
cross section of SN3. Measurements were made at a 500-g load 
for 15-sec duration. Twenty measurements were made at the high- 
pore-fraction and surrounding regions. The error in the 
microhardness measurement was estimated at lOpercent from 
repeated trials. 


5.0 Results and Discussion 

The following sections show representative images from the 
contact-scanned Si 3 N 4 disks. Reflection coefficient, attenuation 
coefficient, phase velocity, and cross-correlation velocity images 
obtained from contact scanning are designated as /?(/), «(/), 
V^f), and Vp U j se images, respectively. 

5.1 Ultrasonic Contact Scan Images and Properties Versus 
Frequency 

Figures 1 0 to 13 show /?(/), a(/), V p (f), and V pu i se images and 
mean values at 30 to 110 MHz in 10-MHz increments for disk 
SN2 containing six distinct regions of high pore fraction. (Nomi- 
nal center frequency was 50 MHz for the transducer used to scan 
SN2.) At most frequencies, a(/), V p (f ), and V pu | se images clearly 
revealed all six regions of high pore fraction. Additionally, V pu j se 
and V p (f) images revealed pore fraction variation between the 
disk central section and edge regions. All image types were more 
“focused” toward the middle of the frequency band (50 to 
80 MHz) and less focused at the frequency band extremes 
(<50 MHz and >80 MHz). (The latter statement is not applicable 
to Vp U i se images, which are formed using the whole pulse contain- 
ing the broad band of frequencies.) This result is consistent with 
the fact that as the frequency approaches (from the peak fre- 
quency) the -6-dB frequencies (50 percent bandwidth) of the 
front-surface, pulse-echo magnitude spectrum (-40 and 90 MHz), 
the signal-to-noise ratio becomes low enough to make the results 
of questionable validity and is responsible for the erratic behavior 
of the ultrasonic properties at the frequency extremes (e.g., the 
decrease in mean /?si 1 N 4 above 90 MHz). Between 50 and 
80 MHz, the spatial patterns of the images did not change 
significantly. The /?(/) images did not reveal the six distinct 
regions at any frequency. Mean values for R(f) showed an 
approximately linear increase with increasing frequency (from 30 
to 90 MHz) and ranged from 0.65 to 0.83 (fig. 10(j)). Mean values 
for a(f) increased with increasing frequency and ranged from 
about 0.2 to 1.2 nepers/cm (fig. 1 1 (j)). Figure 1 1 (k) shows 
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Figure 10.— Frequency-dependent reflection coefficient images and mean values at 30 to 110 MHz in increments of 10 MHz for disk 
SN2 containing six distinct regions of high pore fraction. (Nominal center frequency was 50 MHz for the transducer used to scan 
SN2.) 
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Figure 10. — Concluded. 
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Figure 1 1 . — Frequency-dependent attenuation coefficient images and mean values at 30 to 1 1 0 MHz in increments of 1 0 MHz for 
disk SN2 containing six distinct regions of high pore fraction. (Nominal center frequency was 50 MHz for the transducer used to 
scan SN2.) 
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(k) Observed changes in mean asi 3 N 4 versus changes in 
frequency from 30 MHz. 

Figure 1 1 . — Concluded. 
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Figure 1 2. — Frequency-dependent phase velocity images and mean values at 30 to 1 1 0 MHz in increments of 1 0 MHz for disk SN2 
containing six distinct regions of high pore fraction. (Nominal center frequency was 50 MHz for the transducer used to scan SN2.) 
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Figure 12. — Concluded. 



Figure 13. — Cross-correlation pulse velocity image for disk SN2 
containing six distinct regions of high pore fraction. (Nominal 
center frequency was 50 MHz for the transducer used to scan 
SN2.) 

changes in mean asj 3 N 4 (Aasi 3 N 4 ) versus changes in frequency A/ 
from 30 MHz as a means for comparison with figure 4 and will be 
discussed in section 5.8. Mean values for V p (f) decreased slightly 
(0.3 percent) with increasing frequency (fig. 1 2(j)) indicating this 
Si 3 N 4 to be slightly dispersive. 

5.2. Comparison of Contact Scan Results With Those From 
Other Techniques 

Figure 14 compares imaging results from ultrasonic C-scans 
and microfocus radiography of disk SN6 with a(f = 70 MHz), 


V p (f = 70 MHz), and V pu | se images. First, consider the conven- 
tional NDE imaging methods. C-scan back-surface echo reflec- 
tion images at 1 0 MHz (fig. 1 4(d)) revealed all six distinct regions 
as dark (low ultrasonic transmission) areas. C-scan through- 
transmission images at 10 MHz (fig. 14(f)) did not reveal the 
regions. C-scan surface-wave images at 35 MHz (fig. 14(e)) 
clearly revealed the three front-surface-connected regions of high 
pore fraction as dark (low ultrasonic transmission) areas. 
Microfocus x-radiography (fig. 14(g)) clearly showed all six 
distinct regions. Additionally, dark rings appeared to surround 
and highlight the porous regions in the microfocus x-ray image. 
Neither C-scan nor microfocus x-ray images discriminated the 
more porous central section from the highly dense edge region. 

Consider the images obtained from contact scans. The V pu \ sc 
images (fig. 14(c)) quantified variation from edge to center 
including the variation between the distinct regions of high pore 
fraction and the surrounding central region (~1 -percent velocity 
difference), the variation between the distinct regions of high pore 
fraction and the dense edge region (~4-percent difference), and 
the variation between the more porous central section and the 
dense edge region (~3-percent difference). (The percentage change 
was obtained by dividing the velocity values by the velocity 
expected for a fully dense Si 3 N 4 sample (i.e., the theoretical 
velocity of 1 . 1 2 cm/^isec) (ref. 14)). The magnitude and direction 
of these variations in the V pu | sc image are similar to those seen for 
volumetric pore fraction variations in disk SN2 (fig. 7(b)). The 
V p (f) images quantified variation over a broad band of frequen- 
cies in a similar manner to that observed in the cross-correlation 


16 













(°) ^puise ' ma 9e; nominal center frequency, 50 MHz. (f) Through-transmission C-scan (10 MHz focused). 



L_i 
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(g) Through-transmission 
microfocus radiograph. 

Figure 14. — Comparison of conventional imaging methods with ultrasonic contact scan results. 
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velocity images. However, in some instances, V p (f) images 
showed more “scatter” (fig. 14(b)) (note darker-colored dots 
randomly located), and the variation was not quite as obvious as 
it was in the V pu | sc images. The scatter could not be filtered 
without losing essential image features. The magnitude of veloc- 
ity variation was usually similar in V pu j sc and V p (f) images; 
however, minimum, maximum, and mean values varied slightly 
depending on the frequency of the V p ( f) images and were slightly 
different from those observed for Vp^ images (e.g., mean 
V p (f = 70 MHz) for SN 6 was 1 .065 cm/fisec as compared with a 
mean V pu j se of 1.077 cm/fiscc). 

In the a(f) images obtained over a broad band of frequencies, 
larger attenuation coefficients were observed at the edges of the 
distinct regions of high pore fraction than in the surrounding 
regions (e.g., at 70 MHz for SN 6 an attenuation coefficient of 
~4.5 nepers/cm was observed at the edges of the distinct regions 
whereas the mean attenuation coefficient was 1.04 nepers/cm). 
The central section appeared to have a uniformly lower attenua- 
tion coefficient at 70 MHz than the surrounding edge area 
(~0.1 neper/cm as compared with 0.4 to 1.0 neper/cm). 

The /?(/) images did not discriminate the high-pore-fraction 
regions from the surrounding regions in disk SN 6 . 

53 Comparison of S 13 N 4 Disks 

Figures 15 to 19 show a(f = 60 MHz and / = 80 MHz), 
V p (f = 60 MHz and/= 80 MHz), and V pu j sc images of the SN2, 
SN3, and SN5 disks. Disk SN2 contained six distinct regions of 
high pore fraction, disk SN3 contained one such distinct region, 
and disk SN5 contained no such regions; the a(/), V p (f), and 
K pu ise images clearly represent the differences in the disks. The 
more such distinct regions, the higher the mean attenuation 
coefficient, and the lower the mean phase and cross-correlation 
velocities. The V pu j se images consistently were the most revealing 
of disk nonuniformity and differences between disks. The R(f) 
images did not differentiate the microstructual differences between 
the disks. 

5.4 Correlation of Ultrasonic Properties With Each Other 
and With Microstructure 

Figures 20(a) to (c) show ultrasonic properties versus posi- 
tion for a section intercepting two of the high-pore-fraction 
regions in the SN2 disk. From first glance at a(f = 60 MHz), 
V p (f = 60 MHz), and V puIsc images of SN2 (figs. 15(a), 17(a), 
and 19(a), respectively), it appears that variations in velocity and 
attenuation coefficient correspond with each other and with the 
pore fraction profile shown in figure7(b). In reality, disk posi- 
tions where changes in velocity with position dV/dP were most 
severe corresponded to areas of highest «(/) (fig. 20 (c)), and 
velocity variations corresponded exactly with pore fraction varia- 
tions (compare fig. 7(b) with fig. 20(a)). Note the points on the 
curves in figure 20 (c) at the 6 - and 31 -mm positions; the attenu- 
ation coefficient achieved maximums where dV/dP was most 
negative (at X ~ 6 mm) and most positive (X ~ 31 mm). These 
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(a) SN2. 


(b) SN3. 



(c) SN5. 

Figure 1 5.— a ( f = 60 MHz) images of SN2, SN3, and SN5. 


results are consistent with those from a previous study (ref. 16) 
of ultrasonic imaging of sintered a-SiC disks. In that study, it 
was found that regions in a(f) images showing high attentuation 
coefficients generally overlapped boundaries seen in the corre- 
sponding Vpyisc image (i.e., a high attenuation coefficient was 
observed where large spatial gradients in velocity and density 
existed). A plausible explanation for these results follows: It is 
believed that refractive scattering of waves occurs at imped- 
ance-mismatched material boundaries (see the distinction between 
the high-pore-fraction and surrounding regions in fig. 8 ). This 
results in an apparent attenuation coefficient being measured as 
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(b) SN3. 



(c) SN5. 

Figure 1 6. — a ( f = 80 MHz) images of SN2, SN3, and SN5. 



(a) SN2. 




(c) SN5. 

Figure W.—V p (f = 60 MHz) images of SN2, SN3, and SN5. 
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(a) SN2. 



(b) SN3. 



(c) SN5. 

Figure 18. — V p (f = 80 MHz) images of SN2, SN3, and SN5. 



(a) SN2 (nominal center frequency, 50 MHz). 

cm/pisec 



(b) SN3 (nominal center frequency, 125 MHz). 

cm/^sec 



(c) SN5 (nominal center frequency, 125 MHz). 
Figure 19. — V pu | Se images of SN2, SN3, and SN5. 
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Figure 20. — Ultrasonic velocity versus position for section inter- 
cepting two high-pore-fraction regions in SN2 disk. 

ultrasonic energy is redirected. The effect on the measured 
attenuation coefficient is difficult to predict, since it is difficult to 
know exactly how the ultrasound is redirected. Thus, such re- 
fractive scattering effects may lead to both high and low values 
for the measured attenuation coefficient. In this study, high 
attenuation coefficients were observed at relatively sharp pore 
fraction boundaries. 


It is also interesting to note that the dependence of a(f) on 
frequency was quite different at a relatively sharp pore fraction 
boundary than in a surrounding region in the SN2 disk (fig. 21). 
Therefore, the slope of a(f) (dc^fydf) may also provide infor- 
mation on spatial differences in microstructure (refs. 34 and 35). 

The curves of pulse velocity and phase velocity are similarly 
shaped, but the reflection coefficient variations (fig. 20(d)) do not 
correspond with variations in the other curves of ultrasonic 
property versus position. 



c 



Frequency, f, MHz 
(b) Away from sharp boundary. 

Figure 21 . — a (f) ± a error properties versus 
frequency for SN2. 
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5.5 Comparison of Ultrasonic Contact Scan Images for 
Scans on Machined and Polished Disk Faces 

Figures 22 to 26 show a(f= 60 MHz) and the percent error in 
cc(f= 60 MHz), V p (f= 70 MHz), V pulse , and R(f= 60 MHz) images 
for the experimental configuration where the ultrasonic trans- 
ducer was scanned first on the machined face and then on the 
polished face (opposite surface) for disks SN6, SN3, SN4, and 
SN5. In most cases, no significant differences in feature differen- 
tiation were noted between the same image types for machined 
and polished face scans (note that each disk was flipped over for 
the polished scan after scanning on the machined surface). Mean 
values of a(f = 60 MHz) for SN6 from machined and polished face 
scans differed by about 25 percent, but this is within the mean 
experimental error (~30 percent) for the a(f = 60 MHz) image 
(fig. 23). The minimum, maximum, and mean values for 
Vp(f= 70 MHz) images of SN3 were nearly identical for machined 
and polished face images. More scatter dots (discussed previ- 
ously) were apparent in the machined face image than in the 
polished face image for V p (f = 70 MHz), and it appears that 
slightly more detail was differentiated in the polished face image. 
However, the latter two results were not always observed for 
machined versus polished V p (f) images and are probably caused 
by small scan-to-scan and within-scan couplant variations. Mini- 
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fa) Machined surface. 



(b) Polished surface. 

Figure 22. — a (f = 60 MHz) images for SN6 with transducer on 
machined and polished surfaces. 



(a) Machined surface. 



(b) Polished surface. 


Figure 23.— Percent error in a [f = 60 MHz) images for SN6 with 

transducer on machined and polished surfaces. 

mum, maximum, and mean values for V pu i se images of SN4 were 
essentially identical for the machined and polished face images. 
Mean and maximum values for R(f = 60 MHz) images of SN5 
were each very similar for machined and polished face images. A 
low reflection coefficient strip was apparent in the polished face 
image but not in the machined face image. The latter result was not 
always observed for machined versus polished R(f) images and 
probably results from couplant thickness variations (fig. 3(c)). 
The results of the machined versus polished face scans were 
surprising considering the expected sensitivity of frequency- 
dependent properties on surface condition and allow the reliable 
application of the contact scan method to materials with machined 
as well as polished surfaces. 

5.6 Dry- Versus Wet-coupled V pu | se Images 

In figure 27, a wet-coupled V pu i sc image of disk SN2 (obtained 
at a nominal transducer center frequency of 50 MHz and a 
transducer diameter of 0. 1 25 in. (0.3 cm)) is compared with a dry- 
coupled Vp U \sc image (obtained at a nominal transducer center 
frequency of 10 MHz and a transducer diameter of 0.25 in. 
(0.6 cm)). Both images clearly reveal and similarly quantify the 
six distinct regions of high pore fraction. This result is important 
for materials that cannot be ultrasonically examined with wet 
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(a) Machined surface. 



(b) Polished surface. 

Figure 24. — Vp (f = 70 MHz) images for SN3 with transducer on 
machined and polished surfaces. 



(a) Machined surface. 



(b) Polished surface. 

Figure 26 . — R ( f = 60 MHz) images for SN5 with transducer on 
machined and polished surfaces. 


cm/jtsec 



(a) Machined surface. 


cm/psec 



(a) Wet coupled (nominal transducer center frequency, 50 MHz; 
transducer diameter, 0.125 in.) 


cm/pisec 



(b) Polished surface. 

Figure 25.— y pulse images for SN4 with transducer on machined 
and polished surfaces. Nominal center frequency, 125 MHz. 


cm/p. sec 



(b) Dry coupled (nominal transducer center frequency, 10 MHz; 
transducer diameter, 0.25 in.) 


Figure 27 .— V pulse images of disk SN2. 
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coupling. The high-pore-fraction regions are most clearly revealed 
in the wet-coupled image, as is consistent with the facts that the 
beam diameter for the wet-coupled configuration was approxi- 
mately half that for the dry-coupled configuation and that the wet- 
coupled configuration had approximately an order of magnitude 
higher signal-to-noise ratio than the dry-coupled configuration. 
The minimum, maximum, and mean values were uniformly 
about 1 percent higher in the dry-coupled image than in the wet- 
coupled image. 

5.7 Longitudinal Versus Shear-Wave V pu | S e Images and 
Poisson’s Ratio Maps 


Figures 28(a) and (b) show ultrasonic longitudinal (nominal 
transducer center frequency, 10 MHz) and shear-wave (nominal 
transducer center frequency, 5 MHz) V pu | se images, respectively, 
of SN2 obtained from dry-coupled scans. (The relatively low 
signal-to-noise ratio for the shear-wave, dry-coupled configura- 
tion resulted in data not being acquired at some scan points near 
the left edge of the disk during the shear-wave scan.) Both images 
clearly reveal and similarly quantify the six distinct regions of 
high pore fraction. Velocity variation was about 4 percent over the 
longitudinal-wave image and about 2 percent over the 

shear-wave Vp U | se image. 

Poisson’s ratio vcan be computed directly from longitudinal 
and shear-wave velocities according to (ref. 13) 


v = 




(19) 


where Vi and Vs are the longitudinal and shear-wave velocities, 
respectively. Predictions of mechanical and thermal behavior of 
components (by finite element analysis, e.g.) traditionally assume 
a constant value for Poisson’s ratio. This assumption can lead to 
reduced accuracy of the analysis if in fact “apparent” Poisson’s 
ratio (ref. 14) varies within a component because of micro- 
structural variation. The use of materials in demanding high- 
performance applications where subtle within-sample mechanical 
of such aberrations (refs. 2 and 3). Figure 28(c) shows the 
apparent Poisson’s ratio image of the SN2 disk constructed from 
the longitudinal and shear-wave velocity images through the use 
of equation (7). The Poisson’s ratio image indicates similar 
variations to those shown in the velocity images and quantifies 
variations in Poisson’s ratio across the disk. Apparent Poisson’s 
ratio varied across the disk from -0.26 to -0.29 (-10 percent 
variation) with a mean value of -0.27. The relative error in the 
Poisson’s ratio image was <1 percent as the longitudinal and 
shear-wave V pu | se images used to form the Poisson’s ratio image 
were reproducible to 0.1 percent. The Poisson’s ratio image can 
in concept be used in a NASTRAN finite element analysis 
(ref. 36) by dividing the disk into a finite element grid, using 
different material cards for finite elements having different 
Poisson’s ratios, and inputting the appropriate value of apparent 
Poisson’s ratio on each material card. 


cm/p.sec 



(a) Pp U ise ima 9 a from dry-coupled, longitudinal wave scan 
(nominal transducer center frequency, 10 MHz). 


cm/p.sec 



(b) V pU | Se image from dry-coupled, shear-wave scan (nominal 
transducer center frequency, 5 MHz). 



(c) "Apparent" Poisson's ratio image of same disk constructed 
from longitudinal and shear-wave velocity images. 


Figure 28— V pulse and Poisson's ratio images of SN2. 


Elastic properties such as Poisson’s ratio are most accurately 
calculated from phase velocities, since the elasticity of a material 
may depend on the loading rate (ref. 23). As previously men- 
tioned, velocity determined from cross-correlation is essentially 
a pulse velocity, as the entire wave train (containing a broad band 
of frequencies) is considered in the calculation. Pulse and phase 
velocity will be essentially equal over the frequency bandwidth 
investigated if dispersion (velocity varying with frequency) due 
to the material itself and sample geometry are negligible. Since 
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essentially no dispersion was apparent for this Si 3 N 4 material, the 
Poisson’s ratio values obtained from cross-correlation velocity 
measurements were deemed to be sufficiently accurate and are 
dose to values obtained in other studies (ref. 37). 

5.8 Experimental Versus Expected Results for Frequency 
and Pore Fraction Dependence 

Figure 29 shows the experimental relations for ultrasonic 
properties versus frequency obtained for SN2 with the ultrasonic 
transducer placed directly at the disk center (-0.03 pore fraction 
averaged through the thickness), at the disk edge (-0.01 pore 
fraction averaged through the thickness), and in the high-pore- 
fraction regions (-0.04 pore fraction averaged through the thick- 
ness) (fig. 7(b)). The following comments pertain to the frequency 
range of approximately 30 to 110 MHz. Below and above this 
frequency range, the low signal-to-noise ratio produced results of 
questionable validity as previously explained (section 5.1). 

Phase velocity V$i 3 N 4 decreased with increasing pore fraction 
at any frequency and decreased slightly with increasing fre- 
quency (fig. 29(a)). These results agree well with those predicted 
(fig. 2(c)). Absolute magnitudes of velocity values and magni- 
tudes of changes in velocity values with changing pore fraction 
correlated with predicted theoretical (fig. 2(a)) and previously 
observed empirical (fig. 2(b)) results. 

The reflection coefficient /? S i 3 N 4 had no predictable depen- 
dence on pore fraction and increased with increasing frequency 
for all pore fractions (fig. 29(b)). Regarding the latter result, 
experimental results agreed with those predicted. The /?si 3 N 4 was 
expected to slightly decrease with increasing pore fraction 
(fig. 3(a)); experimental results showing no pore fraction depen- 
dence for /?si 3 N 4 were surprising, since the regions of high pore 
fraction were readily apparent at disk surfaces (fig. 5). The 
behavior of /?si 3 N 4 versus frequency appears to indicate a couplant 
thickness of -0.5 fum (fig. 3(c)). Comparing the reflection coef- 
ficient for Si 3 N 4 calculated from velocities and densities (/?si 3 N 4 
= 0.46) (ref. 13) with the larger frequency -dependent 7?si 3 N 4 
observed experimentally revealed the necessity of using the 
experimentally observed reflection coefficient to determine accu- 
rate attenuation coefficients. 

Attenuation coefficient a$i 3 N 4 increased with increasing fre- 
quency for all pore fractions (fig. 29(c)). The increases were on the 
order of <2 nepers/cm over the frequency range 30 to 1 10 MHz. 
The asi 3 N 4 versus frequency behavior at distinctly different pore 
fraction regions was not consistently different. A variation was 
expected, since grain size distribution did not appear to vary 
between the high-pore-fraction and surrounding regions. How- 
ever, it is difficult to make conclusive statements regarding the 
effects of pore fraction changes on changes in asi 3 N 4 by using the 
methodology of this study. Apparent attenuation coefficient may 
have been measured and imaged because of scattering at rela- 
tively sharp material boundaries as discussed previously (the 
dependences of a(f) on frequency, and the values for asi 3 N 4 were 
observed to be significantly different between pore fraction 
boundaries and regions away from such boundaries as shown in 



Frequency, f, MHz 
(c) Attenuation coefficient. 


Figure 29.— Experimental relations for ultrasonic properties 
versus frequency obtained for SN2 directly at disk center 
(-0.03 pore fraction averaged through thickness), disk edge 
(-0.01 pore fraction averaged through thickness), and high- 
pore-fraction regions (-0.04 pore fraction averaged through 
thickness). 

fig. 21). Further competing material effects (such as second-phase 
variations) may have also been present, affecting «si 3 N 4 - Expected 
changes in «si 3 N 4 due to scattering from pores (fig. 4(c)) were 
compared with observed changes in mean «si 3 N 4 w ith changes in 
frequency (fig. 1 1 (k)). An average pore concentration n 0 of 
approximately 10 8 /cm 3 was indicated by the observed magnitude 
of change (<2 nepers/cm from 30 to 1 10 MHz) (if grain scattering 
and other contributions to «si 3 N 4 were assumed to be negligible 
and pores were approximately l-/im radius). 

5.9 Correlation of Ultrasonic Contact Scan Results With 
Physical Properties 

Hardness. — The average of 20 Vickers microhardness meas- 
urements in each region yielded similar apparent Vickers hard- 
ness numbers (mean ± standard deviation) of 1126 ± 56 for the 
high-pore-fraction region and 1 103 ± 58 for the surrounding less 


25 


porous regions of disk SN3. Thus, Vickers microhardness meas- 
urements were not able to clearly discriminate the pore fraction 
variation. 

Thermal behavior. — Thermal behavior variations correlated 
with pore fraction variations in disk SN6 as shown in four 
successive time frames separated by 0.05 sec in figure 30. The 
three front-surface-connected regions of high pore fraction, as 
well as the generally less dense central section of disk SN6, were 
highlighted by lighter gray scale (higher thermal energy den- 
sity). Variations were most prominent in frames 1 and 2 and 
began to fade out later in time. The initial magnitude of variation 
between highly porous and surrounding regions was on the 
order of 2 deg C. (None of the the six distinct regions of high 


pore fraction were differentiated in thermal diffusivity images 
when using a similar procedure. Average thermal diffusivity was 
0.23+0.03 cm 2 /sec for disk SN6, which agreed well with those 
previously reported (0.22 to 0.26 cm 2 /sec) for sintered Si 3 N 4 
(ref. 38).) 

6.0 Further Discussion 

The potential of the ultrasonic contact scan technique for 
materials characterization was demonstrated in this study. A 
single transducer in combination with spectral analysis allowed 
material characterization in terms of ultrasonic material proper- 
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(b) 0.1 0 sec. 


(d) 0.20 sec. 


Figure 30. — Thermal wave images of disk SN6. 


(c) 0.1 5 sec. 
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ties over a wide band of frequencies. This method results in time 
savings over scanning using discrete frequency techniques. The 
ultrasonic contact scan technique is applicable to flat machined 
materials in which access to only one side is possible and is most 
useful in a materials development program for quantitative char- 
acterization of sample homogeneity. It is also envisioned that this 
method would be useful in a quality assurance laboratory for 
material coupon analysis. The scan system hardware and soft- 
ware can be upgraded for improvements or additions with the 
following envisioned: First, with the development of wheel 
probes and multielement assemblies for scanning (ref. 39), and 
the ever-increasing power of computers for data acquisition and 
analysis, scan and analysis time can be significantly reduced in the 
future. Second, images of additional ultrasonic properties, such as 
group velocity, peak frequency of the material-modulated ultra- 
sonic echo, and attenuation coefficient slope, can be constructed 
by using the digitized and stored ultrasonic waveform data. These 
properties and images can be used for further analysis and 
verification of microstructural condition. 


7.0 Conclusions 

An ultrasonic scan procedure using the pulse-echo contact 
configuration was employed to obtain maps of pore fraction 
variations in polished and machined sintered silicon nitride 
samples in terms of ultrasonic material properties. Ultrasonic 
velocity (pulse and phase), attenuation coefficient, and reflection 
coefficient images were obtained simultaneously over a broad 
band of frequencies (e.g., 30 to 1 10 MHz) by using spectroscopic 
analysis. Wet and membrane (dry) coupling techniques and 
longitudinal and shear-wave energies were used. The essential 
features of this study were as follows: 

1. Velocity and attenuation coefficient images represented 
average through-thickness microstructure for the sample scanned. 
Ultrasonic velocity (longitudinal and shear wave) images revealed 
and quantified the extent of pore fraction variation in the silicon 
nitride disks. The magnitude and direction of variations in cross- 
correlation velocity images were very similar to those seen for 
pore fraction variations. Velocity images were similar for wet and 


membrane coupling. Attenuation coefficient images revealed 
pore fraction nonuniformity due to the scattering that occurred at 
boundaries between regions of high and low pore fraction. The 
contact scan procedure proved more precise and accurate (i.e., 
more quantitative) than ultrasonic C-scans by characterizing 
inhomogeneity in terms of ultrasonic material properties. Spatial 
patterns of velocity and attenuation coefficient images from a 
contact scan were nearly identical at frequencies within the 50- 
percent bandwidth (50 to 80 MHz) but became unfocused at 
frequency extremes. Mean, minimum, and maximum values of 
reflection coefficient and attenuation coefficient significantly 
varied with frequency throughout the entire bandwidth, but 
values for phase velocity varied only slightly with frequency. 

2. Velocity and attenuation coefficient images were nearly 
identical for machined and polished disks, making the method 
readily applicable to flat machined materials. 

3. Maps of “apparent” Poisson’s ratio constructed from longi- 
tudinal and shear-wave velocities quantified Poisson ratio varia- 
tions across a silicon nitride disk. Adirect application of engineering 
ultrasonics is shown, since the Poisson’s ratio map with its grid of 
values can be input to a finite element problem where ideal 
properties were previously assumed. 

4. Experimental velocity results agreed well with those pre- 
dicted as a function of pore fraction and frequency. Experimen- 
tal reflection coefficient results agreed well with those predicted 
as a function of frequency for couplant thickness of approxi- 
mately 0.5 jj m but did not show the expected dependence on 
pore fraction. Experimental results of attenuation coefficient 
changes were similar to those predicted due to pore scattering as 
a function of change in frequency for a pore concentration of 
10 8 /cm 3 . 

5. Thermal wave images of a disk indicated transient thermal 
behavior variations consistent with observed pore fraction and 
velocity and attenuation coefficient variations. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, February 3, 1993 
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Appendix A 

Errors Associated With Pulse-Echo 
Measurements 

The percent error in the pulse-echo velocity V measurement 
was calculated from (ref. 40) 


%ERRy = 100 




( 20 ) 


where A" and AX are thickness and error in thickness, respectively, 
and rand At are time delay and error in time delay, respectively. 
Figure 31 shows the percent error associated with pulse-echo 
velocity measurements. (Error is also dependent on the error in 
frequency-domain signal processing methods used in the calcu- 
lation of ultrasonic properties.) 

The percent error in the reflection coefficient R measurement 
is calculated from 



Figure 32. — Percent error associated with pulse-echo reflection 
coefficient measurements. 


i +Jn 

%ERR* = 100^ -f— (21) 

R (SNR)(/?) 

where SNR is the signal-to-noise ratio of the first front-surface 
echo without the sample present on the buffer rod. It was meas- 
ured by comparing the peak amplitude of FSi(/) to the peak 
amplitude of the noise. The SNR was -250 in this investigation 


(wet coupling configuration). Figure 32 shows the percent error 
associated with pulse-echo reflection coefficient measurement. 
The error in R was low except at low values of R and SNR. 

The percent error in the attenuation coefficient a measurement 
was calculated from 



Figure 31 .—Percent error associated with pulse-echo velocity 
measurements. 
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1 0.5 
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( 22 ) 


where a a is the error in the attenuation coefficient measurement, 
X is sample thickness, Ox * s error in the thickness measurement 

(3 X 10“ 4 cm in this investigation),/? is the reflection coefficient, 
and SNR is the signal-to-noise ratio. Figure 33 shows the percent 
error associated with pulse -echo attenuation coefficient measure- 
ment. Error in a is low for intermediate values of R and low values 
of a. It increases greatly for R< 0.1 and R > 0.9 and as a increases. 
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Percent error 
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Appendix B 

Fabrication of Si3N4 Samples 

Processing 

The microstructure and properties of S^N^based materials 
depend on the processing methods. In this investigation, seven 
near-theoretical-density disk-shaped samples (labeled 

SN1 to SN7) were manufactured from two batches of starting 
powder according to the following procedures (ref. 41). The 
starting powder composition was 91 .7 g Kemanord Si3N 4 + (1 .9 g 
Si02 + 6.2 g Y2O3) densification additives. Two batches of 
powder were milled in ethanol with S^N 4 media in sintered 
reaction-bonded Si3N 4 solid mills for 100 hr. The milled slurries 
were sieved through a 10-jUm sieve and then combined in a large 
plastic bottle. The powders were allowed to settle for about one 
week. Then, excess ethanol was removed by carefully pouring the 
fluid from the container. The bottle containing the slurry was 
shaken and then small samples of the slurry were placed in 
weighing dishes, weighed, dried, and then weighed again in order 
to determine the solids loading (in volume percent) of the slurry. 
A bottle containing slurry was then placed in a hot water bath in 
order to evaporate some of the excess ethanol. A 15 vol % solids 
slurry was obtained. The bottle containing the slurry was tumbled 
for 16 hr. Forty pieces of Si3N 4 media were added to provide 
mixing, but no further particle size reduction was expected. The 
pH of the slurry was raised to 10 by the addition of ammonium 
hydroxide (NH 4 OH). 

Slurry pressing was used to consolidate the samples. Disk- 
shaped samples were formed under 9-MPa pressure in a die. The 
slurry charge was compacted between filter-paper disks, which 
trapped the solids. Porous stainless steel disks were used above 


the top filter paper and below the bottom filter paper for uniform 
pressure application and expulsion of liquid. Aftercareful drying 
to avoid cracking, the disks were sealed with evacuated thin-wall 
plastic tubing. Cold isostatic pressing under 414-MPa pressure 
was performed to increase green density. 

The disks were placed in a tungsten cup with a loose-fitting lid 
and sintered at 2140 °C for 4 hr under 5-MPa nitrogen pressure. 
Disks 1 to 4 were sintered in one run, and disks 5 to 7 were sintered 
in another run. Sintering was done in a tungsten-mesh-element, 
all-refractory-metal, water-cooled, double-wall furnace. High- 
purity boron nitride setters separated the disks from each other 
and from the tungsten cup. 

Machining and Polishing 

Machining was accomplished in two steps. First, a rough 
machining with coolant was performed in which the thickness of 
the samples was reduced from approximately 0.5 cm to 0.3 cm, 
and the top and bottom surfaces were made flat and parallel within 
±25 fd m by using a 20-/im diamond wheel. Then, the samples 
were machined by using the same size diamond wheel so that top 
and bottom surfaces were made flat and parallel to ±2 n m. An 
average of less than 0.02 cm of material was removed from the 
samples in this step. Five samples then were lapped on one side 
to obtain a fine polish by using successively 45-, 30-, 15-, 9-, 
6-, and 3-jUm diamond slurry. Each lapping took approximately 
2 hr. The samples were cleaned with acetone after polishing. 
Figure 9 shows typical surface profiles for machined and polished 
surfaces of the Si3N 4 disks. Machined disk surfaces were charac- 
terized by less than 1-^im peak-to-peak roughness. No difference 
in surface condition between polished regions at a disk center and 
at a high-pore-fraction area were observed. 
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Appendix C 
Symbols 

V L 

Vp(f) 

^pulse 

longitudinal wave velocity 

frequency-dependent phase velocity 

pulse velocity determined from cross-correlation method 

a 

pore radius; scatter radius 

Ks 

shear-wave velocity 

Bit) 

time-domain ultrasonic echo off sample back surface 
where i = 1, 2 

Vo 

X 

velocity in fully dense material 
sample thickness 

Bj(F) 

Fourier transformed B$) where i = 1, 2 

z 

acoustic impedance 

D 

grain size 

a 

attenuation coefficient 

%ERR, 

FS,(0 

percent error in i where i = velocity, reflection coefficient, 
or attenuation coefficient 

time-domain ultrasonic echo off sample front surface 
where i = 1, 2 

«(/) 

7 

Tv 

frequency-dependent attenuation coefficient 
scattering cross section 
normalized scattering cross section 

FS j(F) 

Fourier transformed FS ,(/) where i= 1,2 

e 

phase angle 

f 

frequency 

A 

wavelength 

k 

wavenumber 

V 

Poisson’s ratio 


number of scatterers per unit volume 

P 

density 

P 

pore fraction; position 

o 

error 

R 

reflection coefficient 

T 

time delay between two ultrasonic pulses 

B(J) 

SNR 

V 

frequency-dependent reflection coefficient 

signal-to-noise ratio 

velocity 


time delay between two ultrasonic pulses determined 
from cross-correlation 
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